Radial excitations of Q-balls, and their D-term 



(N 

o 

(N 
C 

m 



43 

Ph. 
CD 
43 



> 
O 
m 

(N 

O 

(N 



X 



Manuel Mai and Peter Schweitzer 

'Department of Physics Yale University, New Haven, CT 06511-8499, U.S.A. 
Department of Physics, University of Connecticut, Storrs, CT 06269, U.S.A. 

(Dated: June 2012) 

We study the structure of the energy-momentum tensor of radial excitations of Q-balls in scalar 
field theories with U(l) symmetry. The obtained numerical results for the 1 < N < 23 excitations 
allow us to study in detail patterns how the solutions behave with N. We show that although the 
fields 4>(r) and energy-momentum tensor densities exhibit a remarkable degree of complexity, the 
properties of the solutions scale with N with great regularity. This is to best of our knowledge 
the first study of the D-term d\ for excited states, and we demonstrate that it is negative — in 
agreement with results from literature on the di of ground state particles. 
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I. INTRODUCTION 

The energy momentum tensor (EMT) is a central 
quantity in the field theoretical description of particles. 
Its matrix elements [l[ give the mass @], the spin Q, and 
the constant d\ of a particle [|| to which we shall loosely 
refer as the D-term. Though not known experimentally, 
d\ is a particle property as fundamental as mass, spin, 
electric charge or magnetic moment. Its physical mean- 
ing is that it gives unique insights into the distribution 
of internal (in hadrons: strong) forces [f|. 

EMT form factors found little practical applications 
0, until it became clear that they can be accessed b 
means of generalized parton distribution functions 0, 
in hard exclusive reactions such as deeply virtual Comp- 
ton scattering Since that the EMT form factors 

were investigated in theoretical frameworks including chi- 
ral perturbation theory, lattice QCD, or effective chiral 
field theories, see 0, H| and fl3l - [l9j . 

Remarkably, in all theoretical studies d\ of pions, nu- 
cleons, nuclei was found negative. A possible explanation 
of this observation provide chiral soliton models [13, EH , 
which describe the nucleon in the limit of a large number 
of colors N c in QCD (20[. In these models the negative 
sign of di emerges as a natural consequence of the sta- 
bility of the nucleon [13, El] . 

To shed some light on the question whether d\ < 
is a general and mo del- independent feature, in Ref. [2l[ 
the EMT of Q-balls was studied. These non-topological 
solitons appear in theories with global symmetries, and 
it is the appearance of the associated conserved charge(s) 
which plays a crucial role for their existence [22l - |24| . 

Q-balls have numerous applications in astrophysics, 
cosmology, and particle physics [25l - |4^ |. They provide 
an extremely fruitful framework for the purpose of clar- 
ifying the relation d\ and stability arguments. In [2l| 
an extensive study of the EMT structure of ground state 
solutions was presented. In all cases d\ < was found, 
and a rigorous proof was formulated that the D-term of 
Q-balls must be negative. Moreover, it was shown that 
stability is a sufficient but not necessary condition for d\ 



to be negative, because some ground state solutions de- 
scribe absolutely stable, others meta-stable or unstable 
Q-balls, depending on the parameters. The general proof 
applies to all cases and always d\ < [2lj . 

This work is dedicated to the study of the EMT of 
radial excitations of Q-balls in scalar field theories with 
U(l) symmetry. To best of our knowledge, this is the 
first study of the D-term going beyond the description of 
a ground state. Radial excitations of Q-balls were studied 
previously in [36j . where the ground state and the first 
two excited states N =1,2 were found for a fixed value 
of the charge Q. In this work, we will work with a fixed 
value of the angular velocity u in the U(l)-space, and 
study the first 1 < N < 23 excitations. With N = 
denoting ground states, the family of Q-ball solutions 
can hence be classified by specifying (Q,N) as done in 
[36|, or by specifying (ui,N) as chosen in this work. 

Our numerical results reach high in the spectrum of 
radial excitations and give fascinating and detailed in- 
sights in the properties of excited Q-balls. In particular, 
we will see that also excited states have a negative d\. 
The present work extends and completes our study of 
the EMT structure of ground state Q-balls. It is impor- 
tant to remark that we make no effort to describe the 
full spectrum of Q-balls which would include also vibra- 
tional or other excitations 12311 . and we will not consider 



quantum corrections [31 1. 

The lay-out of this work is as follows. In Sec. [TT] we 
will briefly introduce the framework, and review how ra- 
dial excitations of Q-balls emerge [36[. In Sec. lIIII we will 
present the solutions for the ground state and radial ex- 
citations 1 < N < 23 which we were able to find with our 
numerical method, and discuss the charge density and the 
EMT densities. In Sec IIVI we will discuss global proper- 
ties like charge, mass, mean square radii, and the D-term 
and investigate patterns how these properties scale with 
N. Remarkably, among the studied quantities d\ varies 
most strongly with N. Finally, in Sec. [V] we will focus 
on the issue of stability and the sign of the D-term. The 
conclusions will be presented in Sec. I VII and some tech- 
nical questions addressed in Appendices. 
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II. Q-BALLS AND RADIAL EXCITATIONS 

In this Section we briefly review the theory of Q-balls, 
and introduce the indispensable formulae on the EMT. 
We use throughout the notation of [HI, and refer to it 
for more details. We study the relativistic field theory of 
a complex scalar field <&(x) with global U(l) symmetry 

C = l(d M *')(d»*)-V, (1) 

where, for suitable potentials V [IH, Q-balls emerge as 
finite energy solutions of the type x) = cxp(io;i) <fi(r) 
with r = |af| and </>(r) satisfying the equation of motion 



cj>"{r) + - <f/(r) + u?<i> - V'{<P) = 
r 

<f>{0) = <h, 0'(O) = 




C<j) b 

2 



(2) 



with 



We will use the potential V(<j)) = A < 
A = 1.1, B = 2.0, C = 1.0 |2l|, El, and set lu' z = 1.37 
which is among the ideal values for our purposes, see 

App.m 

To demonstrate the existence of ground (and excited) 
Q-ball states, one can identify r — >• t and <fi(r) —> x(t) [23| . 
and interpret (J2j> as the Newtonian equation for a unit 
mass particle moving under the influence of the friction 
-Ffrfc = — f i{t) in an effective potential U e g = hui 2 x 2 — V, 



x(t) = F hic - VU oS {x) . 



(3) 



A ground state solution corresponds to the situation that 
the particle starts at t = from rest at xq — > </>o, and its 
motion terminates in the origin x = after infinite time. 

In this picture radial excitations correspond to the sit- 
uation when the particle is given more potential energy 
such that it overshoots the point x = 0, moves "up-hill" 
in the effective potential till it reaches a point of return, 
and finally comes to rest at the origin. In principle, the 
starting points can be chosen such that the particle will 
overshoot the origin 1, 2, 3, . . . , N times, see Fig. [T] 
This means the corresponding solution <f>(r) has N nodes 
at finite r, and we refer to it as the 7V th radial excitation. 
The ground state correspond to N = 0. 

This picture helps to anticipate several features of the 
excitations. As N increases, the particle has to travel 
longer paths, and do more work against the friction. 
Thus we have to release it "close" to the maximum of 
U e s where the effective potential is nearly flat, see Fig.[T] 
The particle has to "wait" there for a sufficiently long 
time before "sliding" down the potential, such the time- 
dependent friction is adequately decreased to allow the 
particle to complete its trajectory. 

Therefore, as N increases, </>o approaches the position 
of the maximum of U e g, see App. [2] and <p(r) ~ <fi 
remains basically constant over increasingly extended 
plateaus "to wait for the frictional force" to diminish. 



The small-r behavior which follows from (|2l is [21 
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FIG. 1: The effective potential U e s(x) = \u 2 x 2 — V(x) as 
used in this work vs. x (thin line). The particle trajecto- 
ries are indicated for N = 0, 1, 2, 3 (solid lines). For better 
visibility for N = (1, 2, 3) the potentials are displaced by 
(0.1, 0.22,0.36) as compared to TV = 0, and the particle tra- 
jectories are displaced by 0.02 after each turn. 



In this Taylor expansion only even powers of r occur, and 
we checked that the coefficients Ck for k = 6, 8, 10, 12 are 
also proportional to U' cS {4>q) though the expressions be- 
come lengthy. This explains why <p(r) exhibits a plateau. 
After the plateau we expect </>(r) to "oscillate" TV-times, 
before it vanishes at asymptotically large r according to 

M 



(r) -> — exp -r^fJl ^; 



(5) 



With our numerical method described in App.lAlwe were 
able to find solutions for the first N = 23 excited states. 

In the following we will discuss the charge density p c h, 
and the EMT densities, namely energy density, T 00 (r), 
pressure and shear force distributions, p(r) and s(r), 
which are given by [HJ 

w 0(r) 2 , 



Pch{r) 
T 00 (r) 



i w V(r) a + ^(r) a 



v{4>), 



s{r) = <p'{r) 2 , 
P( r ) = o^ 2 ^) 2 



G 



0'M 2 - V(4>) 



(6) 
(7) 
(8) 
(9) 



We also define the conserved charge Q = Jd 3 x p c h(r) 
due to the U(l)-symmetry of the theory (U), the mass 
M = Jd 3 x Too(r), and the constant di which can be 
expressed cquivalently in terms of s(r) andp(r) as follows 



di 



1 



M / d 3 xr 2 s(r) 



M 



d 3 x r 2 p{r) 



(4) 



(10) 

The large-r asymptotics ensures that the integrals 
defining Q, M, d\ are well-defined. 
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III. RESULTS FOR THE DENSITIES 

Fig [2] shows the results for the radial fields <f>(r), the 
ground state N = and radial excitations 1 < N < 23. 
The results confirm the features we predicted in Sec. [TTJ 
For N > the initial values 4>o are numerically within 
10~ 6 close to each other. For N > 2 the solutions show 
plateaus with <j>(r) ~ </>o, followed by regions of "oscil- 
latory behavior" with TV zeros, before the exponential 
decays set in according to (0 . For N > 4 the sizes of the 
plateau regions and oscillatory regions are roughly in a 
constant 1 : 3 ratio. 

The N zeros of the solutions <j>{r) imply a strict shell 
structure for the charge distributions p c h{f) which is 
shown in Fig. G3 The N th excited state consists of an 
inner region of nearly constant charge density for N > 4, 
followed by an outer region with N shells. 

Also the energy densities Too(r) in Fig. [4] exhibit char- 
acteristic shell structures. Although they never vanish at 
finite r, the Too(r) show noticeable minima numerically 



very close to the zeros of p c h(r)- This can be under- 
stood in the particle motion picture as follows. We have 
Zoo(r) = §- t E Mn for r e {R t \ <t>{Ri) = 0, 1 < i < N}, i.e. 
the positions Ri, where the fields and hence also charge 
distributions vanish, correspond in time to the transits of 
the particle through the origin, and Tp (i?i) correspond 
to time-derivatives of the kinetic energies at those times. 
In the absence of frictional forces -Ekin would be exactly 
extremal at the origin. Because of friction the extrema 
of .Ekin are somewhat shifted, but those shifts decrease 
with time (-<-> distance) because -Ff r i C oc j. 

For N > 2 the energy densities show "spikes" at the 
edge of the inner bulk region. For N > 3 also the sub- 
sequent inner shells exhibit characteristic "double-spike" 
structures. The reason for that is the contribution of 



the surface energy [21] . The concepts of surface tension 
and surface energy are well defined for ui — > ui m i n (23| . 
but the associated features are noticeable also away from 
this limit [2l[ . If the inner region and the N shells had 
sharp edges, s(r) would consist of (2A^ + 1) (^-functions 
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FIG. 4: The energy densities T 00 (r) as functions of r for < N < 23. 
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FIG. 5: The shear force distributions s(r) as functions of r for < N < 23. 
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FIG. 6: The pressure distributions p(r) as functions of r for < TV < 23. 
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marking the positions of the respective surfaces. For our 
parameters the system is diffuse, but the "smeared out 5- 
functions" in s(r) can be seen in Fig.[5]though the "gaps" 
between the first shells cannot be clearly resolved. 

Also this can be understood in the particle picture, 
where s(r) — > 2E] l - m (t). The zeros of s(r) coincide with 
the turning points in Fig. [TJ The maxima of s(r) occur 
at the positions where the particle is fastest, which is 
close to the origin of the particle coordinate 1 in Fig. [T] 
The characteristic double peaks emerge because the par- 
ticle is slowed down at the origin by the buckle in U e g . At 
earlier times (inner region) the friction -Ff r i c oc j is notice- 
able making the double peaks less symmetric and hard 
to resolve, see Fig. [5] At later times (outer region) the 
friction is diminished, and the double peaks are nearly 
symmetric. 

Fig. |5] show that the pressure distribution of the A th 
excitation changes the sign (2 A + 1) times. Although 
with increasing A the structures are more and more com- 
plex, the results are numerically stable and satisfy the 
stringent tests discussed in App. [B] In particular, in all 
cases the stability condition is satisfied within numerical 
accuracy, as we will discuss in detail in Sec. IVl 

Figs. [2HS] demonstrate that with increasing A the sys- 
tem becomes larger and exhibits an increasing degree 
of complexity. In spite of the complexity, however, the 
size of the system grows with remarkable regularity, as is 
shown in Fig. [7J This Figure displays for the excitations 
1 < A < 23 the respectively first (Ri) and last (Rn) zero 
of the solutions <fr(r). For N = 1 the two radii coincide. 
We observe that the R\ and Rn increase linearly with 
the order of the excitation. 




5 10 15 20 N 



FIG. 7: The positions of the first (iii) and last (Rn) zero of 
the N th radial excitation as function of TV for 1 < TV < 23. 
The discrete data sets are connected by lines to guide the eye. 



To recall, the origin in the particle coordinate x(t) corresponds 
to the zeros of 0(r). The latter are also the zeros of the charge 
distribution and close to the minima of Too(r), see above, which 
emphasizes that all quantities reflect the same shell structure. 



IV. GLOBAL PROPERTIES 

Above we made three important observations which 
will allow us to make predictions for the A-behavior of 
the global (integrated) properties of Q-balls, namely 

(i) the system exhibits a shell structure, 

(ii) the size of the system grows linearly with N, 

(iii) Pc\i{t) and T 00 (r) inside the Q-balls are effectively 
constant independently of A. 

The shell structure of point (i) is evident from Figs. Er- 
|51 The linear growth of point (ii) is apparent from 
Fig. EJ Point (iii) however requires some explanation. 
Strictly speaking the fields <j>(r) and consequently p c h{f) 
and Too(r) are constant only in the inner region, i.e. in 
about 1/4 of the size of excited Q-balls. However, when 
integrating we effectively "average" over the oscillatory 
behavior of these densities in the outer region. There- 
fore, when speaking about global (integrated) properties 
we may think in terms of effectively constant densities 
inside excited Q-balls which motivates assumption (iii). 

On the basis of these observations we expect the fol- 
lowing A-behavior of the charge Q, mass M, constant d\, 
the surface tension 7, surface energy E snr { , and the mean 
square radii (tq), (r%), (r 2 ) of respectively the charge, 
energy, and shear force distributions: 



Q oc N 3 (11) 

M oc A 3 (12) 

di oc A 8 (13) 

7 oc A (14) 

E sulf oc A 3 (15) 

(r 2 ) 1 ' 2 oc A, i=Q,E,s. (16) 

The surface energy is given by E SUI { = Jd 3 rs(r), while 



(tq) = Jd 3 r r 2 p c ii(r)/Q and (r 2 E ) is defined analogously. 
Finally, the mean square radius of the shear forces is 
(r 2 ) = J^°dr r 2 s(r)/j where 7 = J^°drs(r) denotes the 
surface tension. Surface energy and surface tension are 
well motivated notions in the limit u> — > w m i n [23j ] in 
which Q-balls behave like liquid drops [2l|. But they 
will also be helpful in our context. 

On the basis of the assumptions (ii, iii) we expect the 
charge Q and mass M to be proportional to the "volume" 
which grows like A 3 (even though the solutions are too 
diffuse to make "volume" a well-defined concept). The 
scaling predictions (|16[) for the mean square radii also 
follow straight forwardly from assumption (ii). 

The prediction (fl~5"l) for the surface energy is at first 
glance counter-intuitive. One would expect E SUI { to grow 
with "surface area" oc (volume) 2 / 3 oc A 2 . However, we 
have to take into account the shell structure in point (i). 
The ground state has one surface, and the A th excitation 
with its A shells has in addition to that 2 A surfaces. 
The contributions of individual surfaces do grow like A 2 



G 




FIG. 8: Various Q-ball properties plotted vs. N k with the power k chosen according to the predictions in Eqs. (|11H16|I . The 
shown properties X and the corresponding powers k, written as pairs (X, k) are: (a) charge (Q, 3), (b) mass (M, 3), (c) D-term 
(di,3), (d) surface tension (7, 1), the square roots of the mean square radii of the (e) charge distribution ((rg) 1 ^ 2 , 1), (f) energy 
distribution ((r^) 1 / 2 , 1), (g) shear force distribution ((r 2 ) 1//2 , 1), and (h) the surface energy (_E BU rf,3). The discrete data sets 
are connected by lines to guide the eye. 



as the size of the system grows oc N according to point 
(i). But also the number of surfaces grows oc N, which 
yields (|15[i . Similarly, we expect the surface tension 7 as 
defined in [2l], [23[ to be also proportional to the number 
of surfaces, hence the prediction ([14)) . 

In order to derive the scaling behaviour of d\ we may 
use dimensional arguments. The dimensionality of d\ is 
(mass x size) 2 and with mass oc iV 3 and size oc N we 
obtain the prediction (|13|) . Alternatively we may explore 
the liquid drop limit in which df rop = — M 7 i? 4 where 
R denotes the radius of the drop [f| [2l| . With the scaling 
predictions (|12[ [T4l ITS)) for M, 7, and size of the system 
we are again lead to the prediction (|T3|) . 

Fig.|8]shows the global properties Q, M, di, 7, (rf) for 
i = Q, M, s and E auT { plotted as functions of N k with 
the powers k as predicted in Eqs. ([TTIITd]) . The results 
fully confirm the predictions (|11II16I) . Hardly visible in 
Fig.[8]is that for N = 0, 1, 2 the global properties exhibit 
deviations from the scaling behavior (|ll|[To]) . But for 
N > 2 the numerical results follow Eqs. (fTTHTo]) with 
very good accuracy, see Fig. [5] 

In particular, we observe 7 oc TV as predicted in (1141) . 
The more adequate property characterizing the "surface 
tension" at the "boundary" between Q-matter and vac- 
uum is the rescaled quantity 7resc = 

7/ (27V + 1) which 

Y rese W Ar*/<r*> (b) 

0.6 



1 ■ ' ' ■ 1 1 ' ■ ■ ' 1 

5 10 15 20 N 5 10 15 20 N 

FIG. 9: (a) The true "surface density" 7 rcsc = -y/(2N + 1) 
as function of N. (b) The measure of the diffuseness of the 
system Ar 2 /{r 2 ) as function of N. 



takes into account that the N radial excitation has 
(2N + 1) surfaces. Fig. ^ shows that 7 rcsc is nearly 
independent of N as expected. We stress that 7 rosc is an 
average. Q-matter in excited Q-balls does not possess the 
same "surface tension" everywhere, otherwise the peaks 
in s(r) in Fig. [5] would be all equally high. 

We would like to stress that although qualitatively here 
the liquid drop picture is useful, the concept of a surface 
tension is well justified only in the limit ui — > w m j n where 
the solutions exhibit "sharp edges" [23j]. For ground 
states Ar 2 /(r 2 ) can be used as a measure for the dif- 
fuseness of the system, wher e ( Ar 2 ) 2 = (r 4 ) — (r 2 ) 2 
with (r 4 ) = J °°dr r A s(r)/^f [2l|. If for ground states 
Ar 2 /(r 2 ) <C 1 one has "sharp edges" [2lj |. For our pa- 
rameters Ar 2 /(r 2 ) ~ 0.75 for N = 0, i.e. this condition 
is not convincingly realized; the system is diffuse. If we 
apply this measure also to excitations, we find that they 
are similarly diffuse to the ground state, see Fig. [SJd. 

For uj — > w m i n the s(r) would become proportional 
to the sum of (2N +1) 5- functions with support at the 
positions of the surfaces of the shells [2l| . If we assume 
for simplicity the surfaces equidistant and the coefficients 
of 5-functions equal (this is not accurate, see Sec lIIIl but 
will be irrelevant after we take the limit N — > 00 below) 
we would expect that (r 4 ) oc 

Efc=i +1 k 4 /(2N + 1) while 
(r 2 ) oc k 2 /(2N + 1) and 

Ar 2 2 

lim — rf = — for uj -> uj min . (17) 

This corresponds numerically to 0.894 . . . and is remark- 
ably close to the values observed for N > 2 in Fig. [SJd, 
even though our ui is not close to oj m i n . It would be very 
interesting to test the prediction (|17p for ui closer to w m i n . 
But in such situations radial excitations are difficult to 
find numerically, see App. [S] 

The shell structure can also be studied by looking at 
the charge distribution. Since the p c h(r) vanish at the 
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positions where the fields 4>(r) change sign, this allows 
one to define exactly where a shell starts and where it 
ends. The last shell, of course, has no sharp boundary but 
vanishes exponentially according to ([5])- Let us describe 
briefly how the charge is distributed in the largest exci- 
tation N = 23 our numerical method could handle. The 
inner region carries about 16.7% of the total charge of 
this solution, the first shell 1.64 %, and the second 1.58 % 
which is a global minimum. From here on the percent- 
ages carried by the subsequent shells increase gradually 
until the last shell contains 8.8% of the total charge. 

We did not observe regularities other that with respect 
to individual shells, but we found an interesting pattern 
how the charge is partitioned between the inner region, 
and the shell region. Let us define Qmner as the charge 
contained between < r < R\ where R\ denotes the first 
zero of (j>(r), and let Qshciis denote the charge carried by 
all shells, such that Q = Qi nncr + Qshciis- The interesting 
observation is that as N increases Qinnei/Q — > ^ from 
above, while Q s hoiis/Q — > § from below, see Fig. [TUk . 

Of all global properties studied in this work d\ shows 
the strongest variations with N, as it did for ground 
states when uj was varied pit ]. However, when taking the 
dimensionality of d\ into account, see above, one finds 
that the appropriately scaled constant d\ is bound from 
above and below. In |21j the following inequality was de- 
rived for all solutions of the Q-ball equations of motion 



< - 



di 5 
M^|) < 9 



(18) 



In Fig. HUb we see that the radial excitations satisfy the 
inequality (fT5|) . 

Finally let us mention the interesting relation of d\ to 



the relative "wall width" Ar 2 / (r 2 \ 
can be expressed as 



derived in (2l| which 



(r 2 s ) 




N 



oo 

■"mi* 3 



mm J 
5 



(19) 

Again, although 
in we observe in 



where in the last step we used (|17[) . 
in our calculation u> is not close to u>- 
Fig. [TUb that the numerical results are close to the limit 
derived in (fl"9l). 
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V. STABILITY AND di 

For all solutions we find M < mQ where m = w max de- 
notes the mass of a Q-quantum. For the ground state this 
inequality implies absolute stability. But the radial exci- 
tations can decay. For all our excitations lighter ground 
state configurations exist with the same charge. 

For example, our first excited state of uj 2 = 1.37 has 
Q = 342 and M = 461. The following absolutely stable 
ground state solutions have the same total charge but a 
smaller total mass: 

• one Q-ball of us 2 = 0.51 is 1.61 times lighter, 

• two Q-balls of uj 2 = 0.61 are 1.45 times lighter, 

• three Q-balls of uj 2 ~ 0.68 are 1.35 times lighter, 



• fifteen Q-balls of a; 2 = 1.18 are 1.008 times lighter. 

The latter is the threshold for symmetric configurations, 
and 16 Q-balls with uj 2 = 1.21 would be 0.5% heavier. 
Also asymmetric configurations with lower energy exist. 
E.g., the ground states for uj 2 = 0.516 and uj 2 = 1.37 (i.e. 
the groundstate of our excitation) have the same total 
charge but are 1.55 times lighter than the first excited 
state of uj 2 = 1.37. The still heavier excitations N > 1 
have accordingly more decay modes. 2 In short, all radial 
excitations are unstable. 

Nevertheless, the solutions with TV > of course also 
minimize the energy functional, though they correspond 
to local minima of the action. One way to test this offers 
the stability condition, or "von Laue-condition" (45j . 



oc 

/ 



dr r 2 p(r) = 



(20) 



10 15 20 N 



15 20 N 



which was proven to be satisfied for all finite energy so- 
lutions in the Q"ball system in (2l| . It furthermore was 
shown that for all finite energy solutions the pres sure 
is positive for small r and negative for large r [2l|. In 
Sec. IIHI we have seen that the pressure distribution of the 
7V th radial excitation exhibits this pattern and changes 
sign (2^ + 1) times. It is instructive to look in some more 
detail how excited Q-balls realize the condition ([2"l)|) . 

Fig. [TT] shows r 2 p{r) as function of r for the ground 
state and the radial excitations. In spite of the com- 
plexity of the results the condition (|2"U|) is satisfied 
within numerical accuracy which can be quantified as 
follows. For instance, for the ground state we obtain 



FIG. 10: (a) The relative contributions of the inner region 
(circles) and the shell region (triangles) to the total charge 
as function of N. (b) The constant di in units of M 2 (r%) 
(squares) and M E SUT f(r1) (circles) as function of N. 



2 Here we content ourselves to observe that more stable config- 
urations exist, and are not concerned with the dynamics of the 
possible decays. All numbers quoted for up ^ 1.37 are from [2 jl . 
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FIG. 11: r 2 p(r) as functions of r for < TV < 23. Except for the first column the scales on the r-axis are kept constant for a 
better comparison. The shaded regions above and below the r-axis have equal areas such that J^°drr 2 p(r) = 0. 
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FIG. 12: r 4 p(r)/Ci aat as functions of r/i?i ast for < TV < 23, where Ru B t denotes the position of the last minimum of r 4 p(r) and 
Ciast = |^fastP(^iast)|- With these units the global of all curves occur at r/Ru st = 1 and assume the value r 4 p(r)/Ci as t = — 1 
which makes a comparison easier. r 4 p(r) is the integrand of di. The figure demonstrates how the negative sign of d\ appears. 



| f£°dr r 2 p(r)\/ J^°drr 2 \p(r)\ = 0(1O -8 ) and similarly 
up to N < 4. With increasing N it becomes more diffi- 
cult to maintain this accuracy. For N > 10 the accuracy 
is in the range 0(1CT 5 ) to O(10~ 3 ). 

The regions with positive pressure provide forces di- 
rected towards outside. These repulsive forces are com- 
pensated by negative pressure regions with attractive 
forces directed towards the center. Repulsive and at- 
tractive forces cancel precisely according to ([20"]). 

It is interesting to note that the role of the shells is to 
compensate the repulsive forces from the core. In fact, 



on average the shells contribute attractive forces. 

In Ref. [2l[ it was shown that the pattern how the 
pressure distribution satisfies the condition at once 
implies that the constant d\ must have a negative sign. 
Although the sign of d\ can also be deduced from the 
shear forces [2l|, this indicates a connection between d\ 
and stability. In [2l| only Q-ball ground states were stud- 
ied, for which p(r) changes sign only once. Nevertheless 
the general proof that the stability condition (|20| implies 
e?i < in [2l| was formulated assuming that the pressure 
change the sign an arbitrary odd number of times. This 
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is the situation we encounter for radial excitations, and 
our results illustrate how the stability condition (|2T)|) de- 
termines the sign of d\ . 

Fig. shows r 4 p(r) as functions of r. Clearly, inte- 
grating this function over r yields a negative number, and 
up to a prefactor of 5irM the constant d\, cf. Eq. (fT0|) . 
Our results for the pressure distribution therefore fully 
confirm the general proof of the negativ e sig n of d\ from 
the stability relation ([20)1 formulated in |2lj . 

We remark that in the proof of [2l| also the possibility 
was considered that the p{r) could become zero at some 
point without changing sign. We do not encounter this 
situation for the parameters used in this work. 



VI. CONCLUSIONS 

We presented a study of the energy momentum tensor 
of Q-balls in a scalar field theory with U(l) symmetry. 
While in a previous work we investigated in detail ground 
state solutions for different uj [2l|, in this work radial 
excitations of Q-balls were in the focus of our study. 

In Ref. [3(| the radial excitations N = 1, 2 were stud- 
ied previously for fixed charge Q, in other words the exci- 
tations were classified by specifying the charge and order 
(Q, N). Here we adopted a different classification scheme 
and fixed w, i.e. the excitations are specified by (uj,N). 
We were able to find numerically solutions for the ground 
state N = and 1 < N < 23 excitations. All solutions 
obtained in this work were exact solutions of the equa- 
tions of motion. The numerical results were subject to 
stringent tests to guarantee their correctness. On the ba- 
sis of our results reaching high in the spectrum of radial 
excitations we were able to obtain fascinating insights in 
the structure of these excitations. 

As N grows the systems exhibit increasing degrees of 
complexity. The radial field of the N th excitation changes 
the sign iV-times. At the positions Ri with 1 < i < N 
where 1 this happens the otherwise positive 1 charge distri- 
bution vanishes exactly, and the energy density shows at 
positions very close to the Ri clear minima. In other 
words, the charge is distributed over an inner region 
with nearly constant density surrounded by N shells, and 
Too{r) closely follows this pattern. We observed the inter- 
esting pattern that, as N increases, the constant density 
inner region carries 1/6 of the total charge, while the 
remaining 5/6 are distributed on the shells. 

The energy densities show in addition also character- 
istic spikes at the "edges" of the shells due to the impact 
of the "surface energy." The effects of the "surface ten- 
sion" are reflected with even more clarity in the shear 
force distributions. We have shown that the system is 
diffuse for the parameters considered in this work, and 
discussed in which sense the concepts "surface tension" 
and "surface energy" are nevertheless useful. The highest 
degree of complexity is seen in the pressure distributions 
which change the sign (2iV + 1) times. 



In spite of the complexity of the solutions, the proper- 
ties of the excited Q-balls scale with N with great reg- 
ularity. For instance, the size of the system is propor- 
tional to N, independently whether one uses the zeros 
of the <p(r) or square roots of various mean square radii 
to define it. On the basis of general arguments we were 
able to predict also the scaling of other quantities, for in- 
stance M oc TV 3 or di oc N s which are supported by our 
numerical results. Remarkably, among all quantities we 
studied, the D-term varies most strongly with N. Simi- 
larly d\ was the quantity which varied most strongly in 
the study of ground state solutions as functions of u> [21| . 

One of the consequences of EMT conservation is the 
stability (or von-Laue-) condition [H HU stating that 
J °°dr r 2 p(r) = 0. In [2l| this condition was proven an- 
alytically to be satisfied for any solution of Q-ball equa- 
tions of motion, and in this work we could verify numeri- 
cally that the p(r) of radial excitations with its (27V + 1) 
precisely integrates to zero with very good numerical pre- 
cision. 

The important result is that the D-term is negative 
also for all radial excitations. In all approaches where d\ 
was studied so far, it was found negative. But only the 
Z?-terms of ground states were studied so far, and to best 
of our knowledge this is the first time excited states arc 
shown to have also negative D-terms. 

In [2l[ a rigorous proof was given that for Q-balls 
d\ < follows from the stability condition and Q-ball 
equations of motion. In (2~lj only ground state solutions 
were studied for which p(r) changes sign only once. Nev- 
ertheless the proof had to be formulated assuming that 
p(r) could more generally change the sign any odd num- 
ber of times. The results obtained in this work illustrate 
that this is not a pathological case which has to be taken 
into account for the sake of mathematical rigor. Indeed, 
for excited Q-balls one does encounter such a situation 
in practice. 

In this work we also fully confirm another finding of 
21 1, namely that stability is a sufficient but not neces- 
sary condition for d\ to be negative. In fact, we have 
shown that all radial excitations obtained in this work 
are unstable. They correspond to local but not global 
minima of the action, and can decay into configurations 
of absolutely stable ground states with the same total 
charge but a smaller total mass. 

The works presented here and in pTj clearly demon- 
strate the property e?i < for Q-ball systems and, we 
hope, will inspire rigorous proofs of this property also in 
other systems. Our results also establish d\ as a particle 
property particularly sensitive to variations of parame- 
ters of the system. An interesting question remains: can 
d\ be ever positive in a physical system? 
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Appendix A: Technical details 



For that we made the following tests. We checked that 



We assume u> > without loss of generality. Finite 
energy solutions exist for u> in the range [23j 



2V(</>) 



<" 2 < <ax = V"{4>) 



0=0 



For the potential used in this work 0.2 < lu 2 
ground states are absolutely stable for ui 2 < u>. 



< 2.2. 

2 x 
abs ~ 



(Al) 

The 
1.55 

[211 ]. For uj close to uj ni i n it is numerically challenging 
to handle the ground states, let alone radial excitations. 
In order to have an absolutely stable ground state, and 
maximize the chances find numerous radial excitations it 
is profitable to work close to uj 2 hs ss 1.55. In this sense, 
ui = Vl-37 w 0.94 a; a bs is among the ideal choices. 

As N increases, see Sec. lU it is necessary to release 
the particles close to the maximum of U e e given by 




</>constM = 

(A2) 

where the subscript reminds that (|A2[) corresponds to 
one of the "stationary" solutions 4>{r) = const of ([2]) [2lj . 
which however do not satisfy the boundary condition for 
r — > oo. For our parameters C onst = 1-1045 . . . and the 
radial excitations N > 1 are all within 10 -6 of this value. 



Appendix B: Numerical tests 

In view of the complexity of the solutions, it is impor- 
tant to monitor the numerical quality of the solutions. 



A. the stability condition (f20| is valid, 

B. the equation | s(r) + | s'(r) +p'(r) — is satisfied, 

C. the expressions for d± in (fTO]) yield the same result, 



D. p(0) = 2 J°°dr ^p- is equal to p(0) from ©. 



All these relations can be derived from EMT conservation 
[U, QjJ and provide powerful tests for the numerics [2l[ . 
We find relative numerical accuracies between 0(1O~ 9 ) 
and 0(1O -3 ) depending on N and the kind of test. 

In Sec. [V] we already reported how the stability condi- 
tion, test (A), is satisfied numerically. For (B) we checked 
that (2 s(r) + f s'(r) +j/(r))/(2 \s(r)\ + § \s'(r)\ + \p'(r)\) 
is typically of O(10 3 ) or smaller, for r > and V N. 

Concerning test (C): for instance, for the highest exci- 
tation N = 23 we were able to handle with our numerics, 
we obtain from (fT0|) : d\ = —2.0366 x 10 15 using pressure 
distribution vs. df = —2.0360 x 10 15 from shear forces, 
which corresponds to a relative accuracy of 3 x 10~ 4 . 

Concerning test (D): we obtain e.g. for N = - 23 the 
result p(0) = 0.654652 from Eq. while using the 
above quoted formula yields p(0) = 0.654655, which cor- 
responds to a relative accuracy of 5 x 10~ 5 . 

On the basis of these stringent tests we are confident 
that none of the bumps, peaks, structures in Figs. I2TTT21 
are numerical artifacts, but all details of our numerical 
solutions reflect the true characteristics of the excited 
states. 
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